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The determination of the kinetic temperature structure of the troposphere from ground-based 
measurements of oxygen emission spectra in the microwave region is discussed. Molecular absorp- 
tion properties of oxygen and water vapor are reviewed. A new "inversion" technique is described. 
This technique uses a least squares-iteration solution which is applicable to Fredholm integral equations 
of the first kind. The inversion technique is used to reconstruct two types of tropospheric tempera- 
ture profiles. The effects of certain types of errors in the brightness temperature on the derived 
temperature distribution are computed. 



1. Introduction 

In a recent paper by Meeks and Lilley [1963], the 
microwave spectrum of oxygen in the atmosphere 
was discussed, in particular, it was pointed out that 
measurements of the emission spectrum of oxygen, 
such as could be taken from a satellite, could be used 
to determine significant information about the ''gross" 
temperature structure from 40 km down to the earth's 
surface. In the work reported here, the possibilities 
of determining the ''hne" temperature structure of 
the lower 10 km by ground-based measurements of 
the oxygen emission are discussed. The techniques 
discussed here are applicable only for conditions in 
which there are no clouds in the main beam of a highly 
directional antenna. 

In sections 2 and 3, the relevant properties of mi- 
crowave atmospheric radiation and its measurement 
are reviewed, and calculations of the emission spec- 
trum are presented. 

A new method of inverting the radiation measure- 
ments to determine the temperature profile is pre- 
sented in section 4. This inversion technique utilizes 
a variational iterative method of solving integral equa- 
tions, and has the advantage over conventional in- 
version techniques of stabilizing the system against 
errors both in measurement of the radiation intensity 
and in estimating the transmission along the ray 
trajectory. 

In section 5, the inversion technique is applied to 
various atmospheric profiles and the computations of 
the effects of errors in the brightness temperature on 
the derived temperature distribution are presented. 

2. Background 

The amount of radiation flowing in any radiation 
field is conventionally described by the amount of 
radiant energy, dEp, in a specified frequency interval 



(^', v-\-dv) which is transported across an element of 
area da and in directions confined to an element of 
sohd angle do) during a time dt. The sohd angle is 
centered around the polar angles {6, (/)), measured 
with respect to the normal to the surface. This energy 
is expressed in terms of the intensity, h^ by [Chan- 
drasekhar, 1960] 



dEij = h cos Sdadcxidvdt. 



(1) 



The variation of the intensity through a medium is 
related to the absorptive and emissive properties of 
the medium. Absorption is usually expressed either 
as an absorption coefficient, a^, which is the absorp- 
tion per unit length, or as a mass absorption coefficient, 
/c,., which is the absorption per unit mass of absorbing 
material. The two are simply related: 



ai; = Kpp, 



(2) 



where p is the mass density of the material. 

The emission coefficient, 7V, is defined as the amount 
of radiated energy per unit time per unit frequency 
interval per unit solid angle by an emitter of mass 
dm. During conditions of local thermodynamic 
equilibrium, the emission coefficient is related to the 
Planck function, B^Tl by KirchhofPs law [Goody, 1964| 



jv = KvBi;{T), 



where 



BAT) = 



2hv^ 1 

£.2 ghvlkT_l 



(3a) 



(3b) 
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In (3b), h is Planck's constant, k is the Boltzmann 
constant, c is the velocity of light, and T is the abso- 
lute temperature. 

The assumption of local thermodynamic equilib- 
rium (LTE) requires that the distribution of molecular 
energy levels be governed by collisions and that the 
importance of radiative transitions is negligible in 
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maintaining this balance. This imphes that the mean 
time between collisions is much less than the mean 
radiative hfetimes. For the microwave rotational 
transitions of water vapor and oxygen considered in 
this paper, the mean time between collisions is of the 
order 10~^^ sec, while radiative lifetimes are of order 
lO^'^ to 10^^ sec. Thus, at microwave frequencies, 
assumptions of LTE are quite good up to very high 
altitudes. This is to be contrasted with some of the 
infrared and visible absorption bands where assump- 
tions of LTE are valid only in the lowest portions of 
the atmosphere. The reason for this large difference 
between different frequency regions is the cubic de- 
pendence on frequency of the transition probability 
[Schiff, 1955]. 

The expression (3b) for the emission function can 
be simphfied for microwave frequencies by using the 
Rayleigh-Jeans approximation: 



B,{T) = 2kT- 



(4) 



The Rayleigh-Jeans approximation is derived from the 
Planck function by expanding the exponential term in 
the denominator of (3b) in a MacLaurin expansion and 
retaining only the first two terms. The relative error 
in using (4) is about 0.1 percent at 60 Gc/s and 300 °K. 
The amount of radiation emitted by the atmosphere 
and received by a radiometer can be determined by 
integrating the transfer equation along a ray path. 
This equation is [Goody, 1964] 



ds 



-—cxA_h-BAT)l 



(5) 



Here, ds is an increment of path length. The solution 
of (5), evaluated at the ground is 



1 1 -f' 



As')ds' 



Jo 



Bp{s)av{s)e 



Jo 



ds, (6) 



where h,o is the unattenuated intensity of discrete 
sources lying outside the earth's atmosphere. The 
first term represents contributions to the received 
energy from external noise sources attenuated by the 
intervening medium. The second term may be un- 
derstood physically as follows: an element along the 
path emits energy towards earth at a rate equal to 
Bi^aMs; this emitted energy is in turn attenuated by 
the intervening medium by an amount 



exp 



-P aUs')ds' 



the total contribution is then obtained by integrating 
all contributions along the ray path. It is desirable 
to explain the functional dependencies of the vari- 
ables in (6). The absorption coefficient, a^, is an 
explicit (see succeeding section) function of meteor- 
ological parameters such as temperature, pressure. 



and water vapor density. These parameters are, in 
turn, functions of position along the ray path, thus 
making the absorption an implicit function of position. 
It is this dependence only which has been indicated 
above. Similarly, the Planck function, which is an 
explicit function of temperature alone, can be written 
as a function of position. The intensity received at 
the ground, Ip, thus depends on the meteorological 
profile at every point on the ray path. This depend- 
ence is described mathematically by stating that the 
intensity is a functional (as opposed to a function) 
of the profile. The intensity may be related to a 
quantity called the brightness temperature, Thiv): 



h = 2k-Uv), 



(7) 



which compares the received emissions in a single 
direction to those of a blackbody at temperature T^. 

A nonideal antenna placed in a radiating medium 
will receive radiation from all directions for any par- 
ticular antenna orientation, whereas the brightness 
temperature as defined above refers to the contribu- 
tions from only a single direction. The effective 
brightness temperature, Tbi;(0\ 0'), measured by an 
antenna pointing in the direction (0' , (/)') may be ex- 
pressed as [Forward and Rickey, 1960] 



TbAO\ ct)') = -^( G{e\ c/)'; e, ct>)TUO. (t>)dn. 



(8) 



where G{d\ (/>'; 6, 0) is the directional gain function 
in the direction (0, (f)) with respect to the antenna axis, 
Tbu(6, (f)) is the brightness temperature in this direction, 
and the integration extends over all sohd angles, dil. 
For passive probing purposes, a highly directional 
antenna with low side lobes is necessary so that the 
measured effective noise temperature closely approxi- 
mates the brightness temperature of the desired 
direction. This may be expressed mathematically as 
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G(e', (/)'; d, (t>)n,{e, (\>)d^ « Tue\ (/>'), 



where the integration extends over all angles except 
the solid angle of the main beam. Although it has not 
been indicated explicitly, the antenna gain function 
is also a function of frequency, thus further compli- 
cating an already complicated situation. 

The effective noise temperature of an antenna can 
be measured directly by a Dicke radiometer. The 
Dicke radiometer accomplishes this by switching a 
radio frequency signal between the antenna and a 
reference noise source at a constant rate. If a dif- 
ference in noise levels exists between the antenna and 
reference terminals, a "signal" at the switching rate 
is introduced. By knowing the reference source 
temperature and through calibration, the magnitude 
of the difference between the antenna and a reference 
source temperature is found and the antenna tempera- 
ture can be calculated. 
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This minimum detectable signal can be sh own 
[Drake and Ewen, 1958] to be proportional to l/V^r, 
where B is the predetection bandwidth of the receiver 
and T is its integration time. Thus, for measurements 
of high sensitivity, it is necessary to have both a large 
bandwidth and a long integration time. To study the 
atmosphere by radiation measurements, a narrow 
bandwidth and a relatively short integration time are 
desirable. The narrow bandwidth is desirable be- 
cause assumptions of monochromaticity considerably 
simplify calcuations of radiation properties; the inte- 
gration time must be small enough that the gross tem- 
perature and pressure profiles do not change during 
times of measurement. The integration time could 
be increased for periods of relative calm. A band- 
width of about 20 Mc/s and an integration time of 
about 10 sec appear suitable for significant measure- 
ments. Within a frequency interval of this magnitude, 
the source function and the absorption profile remain 
essentially constant and assumptions of monochro- 
maticity are valid. 

3. Gaseous Radiation in the Troposphere 

In the part of the microwave region considered here 
(10 to 100 Gc/s) the principal atmospheric gaseous 
absorbers and, hence, emitters of radiation, are oxy- 
gen and uncondensed water vapor. Water vapor 
absorption is due to a pure rotational transition at 
22.2 Gc/s and the nonresonant contribution of other 
rotational hues starting around 180 Gc/s and extending 
well into the infrared [Van Vleck, 1947b]. Oxygen 
absorption is due to a band of rotational transitions 
starting at about 53 Gc/s and extending to about 66 
Gc/s [Zimmerer and Mizushima, 1961]. Figure 1 
shows the calculated wet, dry, and total absorption 
as a function of frequency. The calculations for a 
standard atmospheric profile were based on the Van 
Vleck equations given below. It is evident from this 
figure that oxygen absorption dominates over that of 
water, except in the immediate region of the water 
vapor line. This occurs because of the large "density 
ratio of oxygen to water vapor, since the mass absorp- 
tion coefficients of the two gases are of the same order 
of magnitude. For frequencies close to the oxygen 
complex, water vapor absorption is much less than 
oxygen, except perhaps for extreme chmatic profiles. 

Water vapor absorption coefficients in the micro- 
wave region may be calculated from the formulas of 
Van Vleck [1947b], 
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Figure 1. Microwave absorption coefficients at standard surface 
conditions. 



In (9), y"20) jg ^i^g absorption (dB/km) at frequency 
p (c/s) due to uncondensed water vapor of density 
p (g/m'^) at temperature T{°K), Vo is the resonant fre- 
quency of the water vapor 6-5 — 5_i electric dipole 
rotational transition (22.235 Gc/s), c is the speed of 
light, and A^' is the hne width, i.e., the half-width of 
the line at half intensity. The line width is equal to 
the reciprocal of the mean time between collisions of 
the absorbing molecule and depends on the pressure, 
temperature, and presence of foreign gases. The 
measured value of the line width at one atmosphere 
pressure and 318 °K is A^o = 2.6 Gc/s [Becker and 
Autler, 1946]. The pressure and temperature depend- 
ence of the water vapor line widths have been given 
as [Barrett and Chung, 1962] 



Ai^ = Aj^o 



1023.25 



318.0V^ 



-V (1-h 0.0046 p). (10) 



In the above expression, the (juantity X as calculated 
by Benedict and Kaplan [1964] for H^O-HjO collisions 
was found to be 0.614. Calculations l)y Benedict and 
Kaplan [1959] for the case of H2O-N2 broadening re- 
sulted in A'n^ = 0.626 between the temperature ranges of 
220 to 300 °K. These calculations indicated a value 
of Aro = 0.09019 cm-i (=2.704 Gc/s) at r=300 °K. 
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Recent calculations [Benedict and Kaplan, 1964] have 
indicated that Xq^ is approximately equal to X^^. 
These considerations indicate, then, that the tempera- 
ture and pressure dependence of the water vapor 
line widths in air may be represented as 



Al^H: 



20 



H20 



""^H^Q I PH2O-O2P02 



p 

PH20-N2^N2 



-^02 



©"'"■ 



(11) 



where Pm is the partial pressure of the mth molecular 
species, y8H20-m represents the relative effectiveness 
of the mlh constituent over H2O — H2O collisions in 
producing water vapor broadening, Tq is some refer- 
ence temperature, and /^(= Ph20 + ^02 ~l" ^N2) is the 
total pressure of the air. As of yet, there has been no 
reason to doubt the linear dependence of the H2O 
line widths on pressure for the range of pressures nor- 
mally occurring in the atmosphere. 

Oxygen absorption in the microwave region is due 
to a band of magnetic dipole rotational transitions 
centered around 60 Gc/s (0.5 cm). The absorption 
coefficient of oxygen is given by [Van Vleck, 1947a] 



CP 

^(02) = ^ ^2 2/^^)^ 



kT 



(12) 



where 



h=lA^ [; 



^vn 



+- 



I^vn 1 



{vN—vf + i!s.v\ {vn-\- 



In (12), C is a constant which depends on the units 
used [Meeks and Lilley, 1963], En is the energy of 
the Mh state, k is the Boltzmann constant, ijln is the 
magnetic dipole moment of the Mh state, vn is the 
resonant frequency of the Mh state, A^ is here a ge- 
neric label for a particular type of rotational transition, 
and the other symbols are as defined above. 

The oxygen line widths, Ai^at, exhibit a more com- 
plicated pressure dependence than do those of water 
[Abbott, 1964], since the ratio ^v IP is not a constant 
with respect to pressure. This behavior is such as to 
make Av/P lower at higher pressures. This has the 
effect of decreasing the nonresonant absorption while 
increasing the absorption at a resonant frequency. 
Measurements have indicated a range of temperature 
dependence of the line widths from about X = 0.1 to 
Z = 0.9 (Ap-T-^^ [Tinkham and Strandberg, 1955; 
Artman, 1953], and uncertainties in line width con- 
stants of about 10 percent. The ratio of O2 — N2 broad- 
ening to O2 self broadening has been measured to be 
0.75 with an uncertainty of 10 percent. 

Figure 2 shows the result of different AvjP constants 
(high- and low-pressure extremes) on the computed 
noise temperature of the atmosphere. Both curves 
were calculated using a linear pressure assumption 
for the line width; the curve corresponding to the 
actual situation will lie somewhere between these two 
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Figure 2. Thermal noise versus frequency for high (0.8 Mc/sjmm 
Ug) and low (2.0 Mc/s/mm Hg) pressure values of the line width. 
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Figure 3. Brightness temperature versus frequency for a vertical 
antenna. 

extremes. 2 These figures indicate that information 
on line widths could possibly be obtained from radia- 
tion measurements of the atmosphere in conjunction 
with simultaneous meteorological measurements 
[Abbott, 1964]. 

A microwave antenna placed in this radiating en- 
vironment will receive energy at a rate which depends 
on the frequency, the antenna orientation, and the 
temperature and density distributions of the atmos- 
phere. For frequencies in the vicinity of the oxygen 
complex, and for low antenna elevation angles, the 



^ All succeeding calculations reported here will be based ( 
ments of Artman [1953]. 



I the low-pressure measure- 
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total integrated absorption (optical depth) is so large 
that the atmosphere is essentially a blackbody radiat- 
ing at a temperature nearly equal to that at the surface 
of the earth. For frequencies near the center of the 
oxygen band, the atmosphere is again a blackbody 
radiating at nearly the surface temperature for all 
elevation angles. In the frequency region from 50 to 
55 Gc/s and for a vertical antenna, however, a large 
change in brightness temperature with frequency 
occurs, as seen in figure 3. Over this relatively small 
frequency interval, a brightness temperature change 
of about 150 °K takes place and in the steepest portion 
of the curve, the slope is about 60 °K per Gc/s. The 
calculations in this figure were based on radiosonde 
data from Dakar, Senegal, February 2, 1956. It is to 
be noticed in this figure that the brightness tempera- 
ture is greater than the surface temperature for fre- 
quencies greater than about 55 Gc/s. This occurs 
because of a temperature inversion which reached a 
maximum of about 294 °K at 1 km. In general, the 
brightness temperature can never be greater than the 
maximum temperature of the atmosphere. 

The possibility is immediately suggested that this 
portion of the spectrum might be useful for a ground- 
based passive probe of the lower atmosphere. On a 
cloudless day, the emission in this region is essen- 
tially a function of the temperature alone because of 
the relatively small water vapor contribution and, 
hence, measurements of this emission might be used 
to infer something about the temperature structure 
of the atmosphere. The large changes in brightness 
temperature, T^, over this frequency interval satisfy 
one of the necessary requirements for a ground-based 
remote probing system; i.e., that the changes in Tb 
with frequency are sufficiently large that a nondegen- 
erate set of equations may be obtained for inverting 
the brightness temperature versus frequency curve 
to determine the temperature-height profile of the 
lower atmosphere. The equipment requirements nec- 
essary to make relevant measurements have been 
discussed in section 2. 

Another important effect takes place in this fre- 
quency region. If the height above which the bright- 
ness temperature contribution may be neglected ^ is 
plotted against the surface absorption coefficient, as 
is done in figure 4, the layers of the atmosphere which 
contribute most strongly to the brightness tempera- 
ture at a given frequency can be determined. In 
this figure it is evident that the relative contribution 
from various altitudes is strongly frequency dependent 
between 52.5 and 55 Gc/s. For frequencies below 
about 52 Gc/s, the radiation contribution is averaged 
over a fairly large height interval (—17 km); above 
55 Gc/s through the center of the oxygen band, the 
contributions to the brightness temperature are from 
a narrow height interval close to the ground. 

The large slope of the brightness temperature versus 
frequency curve and the pronounced changes of the 
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Figure 4. Height above which thermal noise may be neglected 
(99 percent level) versus surface absorption coefficient. 

radiation cutoff height, as pointed out above, arise 
because of the large variations in both the absorption 
coefficient and the optical depth. Thus, in the rela- 
tively small frequency interval, both the low and high 
attenuation limits of the radiation transfer equation 
occur. 

It should be pointed out that both the radiation cutoff 
height and, as has been seen, the brightness tempera- 
ture, depend on the manner in which the fine width 
parameters vary with temperature and pressure. 
Although a more realistic description of this behavior 
will change the numerical values of the quantities 
involved, the basic conclusions of this section will 
remain the same. 

4. Basic Equations of the Inversion Process 

The brightness temperature, r^(j^), at a microwave 
frequency v, observed looking up into the atmosphere, 
is given by [Fryberger and Uretz, 1961] 

Tf,{p) = j^ T{h)a.{h, T{h))e-j/-'^^''-'rO''^^'^^''dh, (13) 

where T is the absolute temperature (kinetic), a:,, is 
the gaseous absorption per unit length, and h is the 
height.^ The problem to be solved is the determina- 
tion of T{h) from a set of measurements of Tbiv). This 
solution is referred to as the inversion process. The 
absorption coefficient a^ is due to water vapor and 
oxygen. As is seen in figure 1, in the frequency re- 
gion 50 to 55 Gc/s, the principal portion of the absorp- 
tion is due to oxygen, and we will hence neglect the 
water vapor contribution in the temperature inversion 
process. It is reahzed, however, that this assumption 
is probably not vahd for very humid conditions, and 
future work should take this into account."^ 

Before deriving the inversion equations, a brief 
review of the relevant concepts and terminology of 
integral equations will be necessary. A Fredholm 



^ This level has here been chosen as the 99-percent level. A more meaningful choice 
could be made if the exact experimental capabilities of the radiometer in use were known. 
For example, if the uncertainty in measurement of Tb were 0.5 °K, then the height level 
corresponding to this value would be chosen. 



* Previously, h denoted Planck's constant. In all succeeding parts of this work, h will 
refer to height. 

^ It is plausible that in conjunction with the oxygen emission measurements, a determina- 
tion of the tropospheric water vapor distribution could be obtained from radiation measure- 
ments taken around the 22.2 Gc/s water vapor Hne. 
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integral equation of the first kind is given by [Tricomi, 
1957] 

(t){x)= J'ix, y)^{y)dy, (14) 

Ja 

where (f){x) arnd^ix, y) are known and ^(y) is the un- 
known function. The function >^'^(;c, y) is called the 
kernel of the integral equation. Equations of the type 
(14) are Hnear, i.e., if ^i is the solution corresponding 
to 01, and ^2 is the solution corresponding to 02, then 
a'^i + fe^2 (a, b are any constants) is the solution of 
a0i H- 602. There are well-known methods for solving 
or getting approximate solutions of this type of linear 
equation if the analytic (exact) form of and.# are 
known [Morse and Feshback, 1953]. Methods of 
solution of nonlinear integral equations are, however, 
not so well known. 

In many cases of physical application, is repre- 
sented by a set of measured values, and no analytic 
form is known for the function. These measured 
values will never be exact because of errors inherent 
in any physical measurement. Such errors in give 
rise to difficulty in the numerical solution of (14) be- 
cause it is unstable, i.e., a small (infinitesimal) change 
in can give rise to a large (finite) change in the so- 
tion '^. This can be shown by the following elegant 
argument given by Phillips [1962]. Consider the 
well-known theorem true for any integrable function 
//(x,y) 



lim fb 



Ja 



equation, and then iterating, correcting the kernel at 
each step [Fow, 1964]. The convergence of this itera- 
tion sequence depends, inter alia, on how the difficul- 
ties in (b) can be eliminated. It was found that 
(before suitable smoothing) reduction of the integral 
equation to a set of simultaneous equations (to be 
solved by matrix inversion) can lead to unreasonable 
estimates of the temperature profile if there are errors 
present in T^. By use of the following least-squares 
method, stability against errors, both in Tb and in 
estimation of the kernel, was achieved. 

To solve (13) for T{h), a standard exponential model, 
a^, is first assumed for the absorption decay with 
height. This linearizes the equation in this first- 
order approximation. Next, the infinite upper limit 
of the integral is replaced by the radiation cutoff 
height, H{v). This height can be very accurately esti- 
mated from a standard model. The spectrum of these 
heights over the chosen frequency range determines 
an upper limit to the height intervals which can be 
probed by this method. Equation (13) is now approxi- 
mated by 

Tb{p)=l T{h)Jf{p,h)dh, (17) 

where the kernel, ;X^(^, h), is determined by the absorp- 
tion model, and is a function of height and frequency 
alone. 

The temperature, T{h), can now be determined by 
inverting (17). A polynomial expansion is chosen 
for T{h), 



y) sin kydy=0 



(15a) 



T{h)=^tiUi{h), 



(18) 



e[x 



'- ^y^'hc, y) sin kydy (15b) 

JQ 



where e(x)-^0 as A:— ><». Equation (15) is sometimes 
called Dirichlet's integral theorem. If (15b) is added 
to (14), there results 



Ja 



(f){x) + e{x) = J'ix, yt'i^iy) + sin ky]dy. (16) 

Ja 

If we now interpret e as an error in the measured value 
of 0, it is seen that a small error can add a finite highly 
oscillatory component to the correct solution ^. Such 
an oscillatory behavior could destroy all physical mean- 
ing of the solution. 

To determine the vertical temperature distribution 
from a set of radiation measurements requires a nu- 
merical solution of the integral equation (13). This 
is difficult because (a) the integral equation is non- 
linear, i.e., the kernel is a function of the unknown 
quantity T{h), and (b) this equation is also unstable 
against errors in the measured quantity Tb(v). The 
difficulty of nonlinearity can be overcome by first 
linearizing the equation by using a standard model 
to estimate the kernel, solving the resulting linearized 



where the set of {u} is to be determined and Uiih) 
is some polynomial of degree i. Substituting (18) 
into (17), 



Tb{v)=J^tiK(v,i), 



(19) 



where 



K{v, i) = J^ Ui{h)a['\h)e -[ -T^''^<^^' dh. 



Rather than solve (19) in its unstable form by matrix 
inversion, consider the quantity /, defined as 



7-2 Tb{vj)-^tiK{vj,i)\ 



N+Nj, 



+ 2 2 yki^tiUihjd-TM. (20) 

/c=iV+l M=0 / 

where the 7-summation ranges over frequencies, the 
^-summation ranges over the known boundary condi- 
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tions on the temperature at the height Aa, and jk are 
undetermined Lagrange multiphers. The quantity /, 
except for the constraint terms in the summation over 
k, may be understood as follows. Each term of the 
frequency summation (over j) can be recognized as 
the square of the difference of the "measured" bright- 
ness temperature, Tb(vj) and the computed brightness 
temperature at that frequency. The computed bright- 
ness temperatures are determined completely when 
the set {ti} is given. It is reasonable to assume that 
the set of these coefficients which minimizes I is the 
best approximation to the temperature profile for the 
given set of radiation measurements and choice of 
polynomial degree. Thus, it is required that the 
variation of / with respect to the parameters, {ti}, be a 
minimum. Setting the appropriate partial derivatives 
of / with respect to tk equal to zero yields, for N^ fre- 
quencies, Ni) boundary conditions and an Mb degree 
polynomial. 



2M/A+ ^' yjUk{h,) = Dk 
A-0, 1, . 



(21) 



i=A'+l 



TV 



where 



and 



Aki = Aik = ^K{vj, i)K{vj,k) 



Dk = ^Uvj)K{vj,k). 

Equations (21), together with the Nh constraining equa- 
tions on the temperature, suffice to determine the set 
{ti}. Equations (21) would be modified if the deriva- 
tives of the /;f functions were included. Tliese deriva- 
tives are small, since the kernel is not a sensitive 
function of temperature, and hence are neglected. 

An iteration process can now be set up in the follow- 
ing manner. From the temperature distribution 
derived above and the relation [Hewson and Longley, 
1944], 

Pj \TJ I' 



(22) 



which determines the pressure between any two layers 
in the atmosphere which have a constant lapse rate, 
a new set of absorption coefficients and kernel func- 
tions can be obtained from (12) as a function of height in 
the atmosphere. In (22), M is the molecular weight 
of air, g is the acceleration of gravity, R is the general 
gas constant, and In is the lapse rate between layers 
1 and 2. A numerical integration is then performed 
using the derived set of absorption coefficients to 
determine a new set of functions K{p, i) of (19), taking 
fully into account the temperature and pressure de- 
pendence of the absorption coefficients and the line 
widths. Equations (21) are again used to determine 
the {ti} by utihzing the new estimate of the kernels. 
The iteration process is then continued until a self- 
consistent set of temperature coefficients is obtained. 



By self-consistency, it is here meant that neither the 
kernel functions nor the set {//} changes from the 
Nlh to the (N-\- l)th iteration. 

The numerical integration of the kernel functions 
of (19), for a polynomial of arbitrary degree and with 
exponential absorption within any level, is nontrivial. 
The details of this integration are shown in appendix A. 

It should be mentioned that the least-squares method 
presented above is quite general and could be used 
for other types of assumptions for the temperature 
distribution such as exponential, trigonometric, con- 
stant layer, etc. The number of parameters that 
can be determined accurately is, of course, basically 
limited by the inherent instability of the integral 
equation. 

5. Application of Inversion Method 

In this section, the inversion technique described 
in section 4 is applied to atmospheric temperature 
and pressure profiles from Denver, Colo., and Dakar, 
Senegal. The brightness temperatures were calcu- 
lated using the radiosonde data from the two stations; 
from this basic set of numbers, the inversion method 
of section 4 was apphed in an attempt to reconstruct 
the original profile. The Denver profile was chosen 
because of close similarity to a standard profile; 
Dakar illustrates the effects of a strong temperature 
inversion at low altitudes. 

In figures 5 and 6, the inversion of the August 16, 
1953, Denver, Colo., profile is shown. The inversion 
was accomplished by using 12 frequencies to fit a 
fourth-degree polynomial which was constrained 
below by the surface temperature and above (at 16.0 
km) by a standard atmosphere temperature. The 
derived profile of the lower atmosphere is insensitive 
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Figure 5. Absorption profile as determined by inversion of the 
brightness temperature equation for profile of Denver, Colo. 
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Figure 6. Temperature profile as determined by inversion of 
brightness temperature equation for a Denver, Colo., profile. 

to the choice of this upper constraint. The iteration 
sequence converged to the indicated values on the 
fifth iteration. It is to be noted that convergence 
occurred even though the first approximation, as deter- 
mined by a crude absorption model, is very far from the 
correct values. It is also to be noted that, after a 
sufficient number of iterations, the computed absorp- 
tion profile is almost coincident with the actual, even 
though there is a residual error in the temperature 
prediction. This is due to the relative insensitivity 
of the absorption coefficient to small temperature 
changes. 

The errors in the computed temperature structure 
up to the Nth height level may be expressed in terms 
of the rms error in T, S.E. (7), defined by 



S.E.(T)=J-^(T{hi}-nhiW, 



(23) 



where T{hi) and T^{hi) are the actual and computed 
temperatures, respectively, at the iih height level, 
hi. For the profile just discussed, the rms error in 
r up to 11.6 km was 41.8 °K for the first iteration, and 
was reduced to 2.6 °K after five iterations. The cor- 
responding errors for the pressure were 27.6 and 1.6 
mbar, respectively. Table 1 gives the calculations 
of the rms errors in T and P up to the 11.6-km level 
at each step in the iteration process. 

Table 1. The rms error in temperature and pressure up to 11.6 km 
05 a function of iteration number for Denver, Colo., August 16, 1953 





Iteration 


rms error in T 


rms error in P 




°K 


mbar 


1 


41.8 


27.6 
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34.1 


21.0 


3 


23.1 


13.4 


4 


10.2 


5.3 


5 


2.5 


1.6 



The method of this paper is also capable of deter- 
mining temperature inversions, as may be seen in 
figure 7. The Dakar, Senegal, profile of February 2, 
1956, was chosen to illustrate low-altitude oscillatory 
temperature behavior which occurs in some areas. 



t- 5 



4 - 



1 

\ 


1 1 1 1 

\ 


1 1 1 


1 1 1 1 1 


1 




To = 290. 36 °K 




s 


\\ 




Po= 1007.0 mbar 




TRUE 








\ \ 


\ 4 th ITERATION 




1st ITERATION^ 


X - 


- 




V\ ^4th ITERATION 






\V\,lst ITERATION 


- 




TRUE-^A\\ 




1 1 1 1 


4 th ITERATION ^ll). 

/TRUE 

1 1 : 1 1 1 1 1 1/ 1 



220 



230 240 250 260 270 

TEMPERATURE IN °K 



290 



300 



Figure 7. Temperature profile determined by inversion of bright- 
ness temperature equation for Dakar, Senegal. 

For seven frequencies and a fifth-degree polynomial, 
an rms error of 2.2 °K and 0.507 mbar in temperature 
and pressure was obtained up to the 10.4-km level 
after four iterations. 

It is also important to determine the effect of vary- 
ing the degree of the polynomial on the error in T and 
on the rate of convergence. In table 2, the rms error 
in T and P are shown as a function of iteration number 
for a third, fourth, and fifth-degree polynomial. 

Table 2. The rms error in T and P up to h = 10.4 km as a 
function of iteration number and polynomial degree for 
Dakar, Senegal, February 2, 1956. 
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The rates of convergence for this profile are not as 
spectacular as those of the Denver profile because a 
much more realistic first approximation for the absorp- 
tion coefficients was used (an exponential model with 
statistically determined scale heights was the 0th 
approximation). It is easily seen that iteration does 
not improve the temperature solution for a low-degree 
polynomial (when there are not a sufficient number 
of adjustable parameters) over a best statistical esti- 
mate. For higher degree polynomials, a two- or 
three-to-one improvement can be obtained by iteration. 
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There is an upper limit, however, to the polynomial 
degree for a fixed number of frequencies which can 
be successfully used. This limit is imposed by the 
presence of errors both in the brif2;htness temperature 
and in the first estimation of the absori)tion. 

Finite errors will be present in any physical measure- 
ment of brightness temperature. Any such error in 
Tt) will induce errors into the derived temperature 
profile. It is important therefore to have some esti- 
mate of the amount of measurement error in Tb which 
is tolerable for a physically meaningful temperature 
distribution. In the calculations reported here, two 
different types of errors were introduced, alternating 
and constant. The first type was introduced by adding 
{—Y\bTb\ to the n\\i T^ (the T^s were arranged in as- 
cending order). This set is referred to as type a. 
The set obtained by adding {—Y^^\bTi)\ to the nxh Tb 
is referred to as b. In figure 8 are shown the tempera- 
ture distributions computed from the a and b sets of 
\dTt)\ = 1.25 °K. The maximum error in the tempera- 
ture distribution was about 10 °K for this magnitude 
of \dTb\. The rms errors of a and b up to the 10.4- 
km level were 7.6 and 6.6 °K, respectively. Table 3 
shows the results of various error calculations. It is 
to be noted that the iteration diverged for \dTt)\ = 2.5 
°K, case a, and that the effect of errors increased 
rapidly above \8Tt,\ = 1.0 °K. It is also important that 
no essential difference exists between the exact and 
0.5° cases. This indicated that meaningful atmos- 
pheric information can be obtained from radiation 
measurements that possess a sufficiently small, but 
finite, error. 

The effect of constant errors in Tt, on the derived 
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Figure 8. Effect of alternating errors of magnitude 1.25 °K in 
brightness temperature on the reconstruction of the temperature 
profile. 

Case (a)— errors of (— )"|6r6| added to the n\h brightness temperature (arranged in ascend- 
ing order). Case (b) — errors of (—)"-^^\8Tb\ added to the nth brightness temperature. 
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temperature profile is not as serious as ''alternating" 
errors as discussed in the last paragraph. For an 
"alternating" error of magnitude 2 °K, the iteration 
sequence was near its radius of convergence; for a 
constant error of the same magnitude, the rms errors 
up to 10.4 km are only 6.1 and 4.3 °K for positive and 
negative errors, respectively. The derived profiles 
in the presence of these errors are plotted in figure 
9. Similar errors of + and — 1 °K gave rms errors in 
temperature of 4.0 and 2.3 °K, respectively, up to 10.4 
km. A much more meaningful analysis of the effect 
of errors in Tb on the derived temperature structure 
could be achieved by introducing in some suitable 
fashion random errors in the set of brightness tempera- 
tures. This was not attempted here because of exces- 
sive computer cost. 

Table 3. Results of error calculations 
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Figure 9. Effect of constant magnitude errors in the brightness 
temperature on the reconstruction of the temperature profile. 
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6. Conclusion 

The mathematical method of inversion presented 
here appears capable, under certain restrictions, of 
recovering the thermal profile of the troposphere from 
a set of sufficiently accurate microwave (50 to 55 Gc/s) 
radiation measurements of a cloudless atmosphere. 
A tolerable amount of measurement error in the bright- 
ness temperature appears to be of order 0.5 °K. 

The principle hmitations (or approximations) of the 
technique appear to be the neglect of the water vapor 
emission (not vahd under humid conditions) and the 
neglect of scattering (not valid when clouds are present 
in the main beam of a highly directional antenna). 
The contribution of water vapor could possibly be 
accounted for by simultaneous radiation measurements 
on the 22.2 Gc/s water vapor line. The uncertainties 
inherent in line width calculations can be at least 
partially removed by simultaneous measurements of 
the brightness temperature and the atmospheric pro- 
file above the radiometer. 



7. Appendix A 

The quantity of interest, K[v,i), defined by (19), 
may be expressed as linear combination of functions 
of the form 



M{. 



n) = /i" a\ 
Jo 



{h)e 



-j\(h')a 



dh. 



(Al) 



For notational convenience, the subscript v has been 
omitted. We wish to evaluate (Al) under the assump- 
tion that the nonzero a(h) may be interpolated expo- 
nentially within any height interval: 



a{h)=^a{hi-i)e h! ' , hi-i ^ h^ hi, 



(A2) 



^here 





1 




ho = 0. 


and 


hf/ =H 


We write 




M{n) = y 


r^ k"o 



In 



aihi) 



hi — hi- 



ih)e' 



'dh 



N 

1 



= S e-^0-i + «^Vi^^J-''j-i/^j^g-(Ai,7), (A3) 



/here 



Tj = J a{h-,)Hi(\-e '«." ),;•= 1, 2, . . . ,N, 



T,,= 



and 



^n,j) 



= ' h'>e "'. 
■'"j-i 



h-h,-i 



"j-v 



dh, 72 = 0, 1, 



The functions g{n,j) presumably can not be expressed 
in closed form. Instead, we express them as an in- 
finite sum of integrals, each of which can be evaluated 
in closed form. This sum has been found to converge 
rapidly for ^ a: ^ 10, ^ // ^ 20. The integral 
defining g(n, j) is of the form 



/ 



r/3 

J a 



dx. 



(A4) 



By using a MacLaurin expansion for the second ex- 
ponential term, we arrive at 



^=1 



b'Ai,n 



(A5) 



where 



i,n=\ X^'i 



{a-ci)x 



dx. 



The Ai, ri are easily obtainable from any standard table 
of integrals. The series defined by (A5) converges for 
any finite value of b or Ai^ „ as may be seen by compari- 
son with the series for e^-^. 

Thus, by making the correspondence a-^/ij-i, 
13-^ hj, a-^-llHj, h -> a{hj-x) Hje^^j-i^^j, and c^ HHj, 
we have succeeded in evaluating (Al). 
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